AFFDL-TR-68- 136 


ON  THE  INFLUENCE  OF  INITIAL  GEOMETRIC 
IMPERFECTIONS  ON  THE  BUCKLING  AND  POSTBUCKLING 
BEHAVIOR  OF  FIBER-REINFORCED  CYLINDRICAL  SHELLS 
UNDER  UNIFORM  AXIAL  COMPRESSION 


N.  S.  KHOT 


TECHNICAL  REPORT  AFFDL-TR-68-138 


OCTOBER  1968 


D  D  C 


'  B 

This  document  has  been  approved  for  public 
release  and  sale;  its  distribution  is  unlimited. 


AIR  FORCE  FLIGHT  DYNAMICS  LABORATORY 
AIR  FORCE  SYSTEMS  COMMAND 
WRIGHT-PATTERSON  AIR  FORCE  BASE,  OHIO 

Reproduced  by  Iho 

CLEARINGHOUSE 

lot  fedora!  kionhfic  &  Technical 
Information  Springfield  '  a  2  2  ?  51 


NOTICE 


Whan  Government  drawings,  specifications,  or  other  data  are  used  for  any  purpose 
other  than  in  connection  with  a  definitely  related  Government  procurement  operation, 
the  United  States  Government  thereby  incurs  no  responsibility  nor  any  obligation 
whatsoever;  and  the  fact  that  the  Government  may  have  formulated,  furnished,  or  in 
any  way  supplied  foe  said  drawings,  specifications,  or  other  data,  is  not  to  be  regarded 
by  implication  or  otherwise  as  in  any  manner  licensing  the  holder  or  any  other  person 
or  corporation,  or  conveying  any  rights  or  permission  to  manufacture,  use,  or  sell  any 
patented  invention  that  may  in  any  way  be  related  thereto. 

This  document  has  been  approved  for  public  release  and  sale;  its  distribution  is  un¬ 
limited. 


Copies  of  this  report  should  not  be  returned  unless  return  is  required  by  security 
considerations,  contractual  obligations,  or  notice  on  a  specific  document. 


400  -  D*cembs?  1968  -  C0455  -  59-1354 


AFFDL-TR-68-136 


ON  THE  INFLUENCE  OF  INITIAL  GEOMETRIC 
IMPERFECTIONS  ON  THE  BUCKLING  AND  POSTBUCKLING 
BEHAVIOR  OF  FIBER-REINFORCED  CYLINDRICAL  SHELLS 
UNDER  UNIFORM  AXIAL  COMPRESSION 


N.  S.  KHOT 


This  document  has  been  approved  for  public 
release  and  sale;  its  distribution  is  unlimited. 


AFFDL-TR-68-136 


FOREWORD 

This  report  Is  prepared  as  a  part  of  in-house  effort  under  Project  No.  1473, 
"Exploratory  Development  in  Structural  Mechanics, "  Task  No.  147306,  'Stress 
and  Stability  Analysis  of  Hetergeneous  Anisotropic  Plates  and  Shells  and  Arches." 
The  work  was  carried  out  in  the  Advanced  Theory  Group  and  Structural  Synthesis 
Group  of  the  Theoretical  Mechanics  Branch,  Structures  Division  of  the  Air  Force 
Flight  Dynamics  Laboratory,  Air  Force  Systems  Command,  Wright-Patterson 
Air  Force  Base,  Ohio.  Dr.  N.  S.  Khot  (FDTR)  was  the  Project  Engineer. 

This  report  covers  work  conducted  during  the  period  September  1967  tc 
May  1968.  It  was  released  by  the  author  in  August  1968. 

This  technical  report  has 
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ABSTRACT 

The  effect  of  Initial1  geometric  imperfections  on  the  buckling  and  £w«$>uckliug 
behavior  of  composite  cylindrical  shells  are  subjected  to  uniform  axial  compression 
is  studied  in  this  report.  The  solution  is  obtained  by  employing  von  Karman-Donnell 
nonlinear  strain-displacement  relations  and  the  principle  of  stationary  potential 
energy.  Numerical  results  are  given  for  various  fiber  orientations  in  the  three- 
layer  shall  consisting  of  either  glass-epoxy  or  boron- epoxy  composites,  with 
different  initial  imperfections.  Results  indicate  that  the  boron-epoxy  composite 
shells  are  less  imperfection  sensitive  than  the  glass-epoxy  composite  shells. 
Isotropic  shells  are  found  to  be  more  Imperfection  sensitive  than  composite 
shells.  It  is  noticed  that  the  increase  or  decrease  in  the  classical  buckling  load 
with  change  in  fiber  orientation  is  generally  accompanied  by  a  decrease  or  In¬ 
crease  in  imperfection  sensitivity  of  the  shell. 
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SECTION  I 
INTRODUCTION 

This  report  is  concerned  with  the  influence  of  initial  geometric  imperfections 
on  the  buckling  and  postbuckling  behavior  of  laminated  anisotropic  cylindrical 
shells.  It  has  been  established  that  geometric  imperfections  reduce  the  buckling 
strength  of  isotropic  cylindrical  shells  and  the  extent  of  the  reduction  has  been 
studied  by  many  authors  (References  1  through  8).  A  similar  reduction  in  buckling 
strength  due  to  initial  geometric  imperfections  is  expected  in  the  case  of  the 
anisotropic  cylindrical  shells.  The  increasing  use  of  fiber  composite  cylinders 
in  aerospace  structures  prompted  this  study. 

von  Karman-Donnell  large-displacement  relations,  modified  to  include  geo¬ 
metric  imperfections,  are  used  in  the  analysis.  The  solution  to  the  problem  is 
obtained  by  the  application  of  the  principle  of  stationary  potential  energy. 

The  importance  of  initial  geometric  deviations  on  the  buckling  behavior  of 
an  isotropic  shell  has  been  recognized  for  some  time.  This  can  be  seen  from 
the  early  investigations  of  Flugge  (Reference  1)  in  1932  and  Donnell  (Reference  2) 
in  1934. 

However,  Flugge's  linear  and  Donnell's  nonlinear  analyses  failed  to  explain 
the  magnitude  of  difference  between  the  experimental  and  the  theoretical  results. 

In  1950,  Donnell  and  Wan  (Reference  3)  changed  the  procedure  adopted  by  Donnell 
sixteen  years  earlier  and  proposed  a  new  method.  This  method  was  followed  by 
several  investigators  (References  4  through  7)  with  some  modifications. 

In  the  procedure  developed  by  Donnell  and  Wan,  and  followed  by  others,  the 
total  potential  energy  of  the  system  is  minimized  with  respect  to  the  parameters 
defining  the  shape  of  initial  deviation  of  shell  geometry.  This  implies  that  the 
shape  of  the  cylinder  is  altered  during  minimization  of  the  total  potential  energy. 
Madsen  and  Hoff  in  Reference  8  were  the  first  to  point  out  the  error  in  this  approach. 
In  this  investigation,  minimization  is  carried  out  only  with  respect  to  terms  of 
the  additional  displacement  in  radial  direction  caused  by  the  applied  load. 
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The  general  theory  of  anisotropic  shells  was  developed  by  Ambartsumyan 
(Reference  9)  and  Dong  et  al.  (Reference  10).  Tasi  et  al.  (Reference  11)  and 
Holston,  et  al.  (Reference  12)  investigated  the  buckling  strength  of  filament  wound 
cylinders  where  they  used  the  small-deflection  analysis  presented  by  Chang  and 
Ho  (References  13  and  14).  The  effect  of  heterogeneity  on  the  stability  of  com¬ 
posite  shells  subjected  to  axial  compression  is  investigated  by  Tasi  (Reference  15). 
Thurston  (Reference  16)  outlined  a  large-displacement  analysis  of  filament  wound 
cylinders  under  axial  compression,  where  he  proposes  to  solve  the  two  non¬ 
linear  partial  differential  equations  by  Newton’s  method.  The  effect  of  fiber 
orientation  on  budding  and  postbuckllng  behavior  of  geometrically  perfect  fiber- 
reinforced  cylinders  under  axial  compression  is  investigated  by  Khot  (Reference  17). 

The  numerical  results  presented  in  this  report  are  for  the  geometrically 
Imperfect  three-layer  glass-epoxy  and  boron-epoxy  composite  cylindrical  shells. 
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SECTION  n 

PROBLEM  FORMULATION 

The  coordinate  system  employed  on  the  reference  surface  is  shown  in  Figure  1; 
x  is  taken  parallel  to  the  generator,  y  is  measured  along  the  circumference,  and 
z  is  positive  in  the  inward  direction. 

1.  STRAIN-DISPLACEMENT  RELATIONS 

The  reference  surface  nonlinear  strain-displacement  relations  based  on  the 

work  of  DonnellJReference  2)  for  a  cylindrical  shell  of  radius  R  with  initial  radial 
imperfections  W  are  given  by 

«.  *  u.«  *  T(“,  >*  -  -F  <»,.  >* 

*»  *  V,Y  +  T(V*  -  ■? (",» >*  -  -S-(W-W)  C) 

*«>'  V,«  +  U,y  +  ",x 

where  U  ,  V  ,  ami  W  are  the  components  of  the  reference  surface  displacement 
vector,  and  W  is  the  initial  deviation  of  the  shell  in  inward  direction  from  a 
circular  cylindrical  shape.  W  is  the  total  displacement  in  the  radial  direction 
including  w  .  Then  ( W  -  W )  is  the  displacement  produced  by  the  applied  load 
only.  The  total  strains  are  then  given  by 

£x  =  ex  "  zkx 


£y  =  €y  ~  zky 


where 


^*y"  exy  “  zk*y 


*«  *  w,««  - 


S  5  w,yy  -»,v 


2(W,«,-V 
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2.  STRESS  RESULTANTS 


The  rels&osis  between  resultant  stresses  and  strains  for  the  entire  laminated 
wall  are  written  as  (References  10  and  18) 


[»]■  MM -MM 
M-  MM  -MM 


where  the  elements  d(3s3)  matrices 


M-M 


are  given  by 


2  c(kl(fe  -h  ) 

k=!  «l  tt+l  k 

-|-2  Cp*(h  *  -hj  ) 
2  b=|  ij  k  +  l  li 


1  y  h  s  .s  . 

■?£,  Cij  ‘Vm*  V 


til 

In  the  above  definitions  n  is  the  number  of  lamina? ,  k  denotes  ihe  k™  layer 
bounded  by  z  -  coordinates  and  hfe4.j  (Figure  1(b)),  and  is  the  aniso¬ 

tropic  stiffness  matrix  associated  with  the  k**1  lamina.  Inversioa  of  Equation  4 
and  sebstitutten  Into  Equation  5  yield 


M-MM  *MM 
[“]-  MM  -  KIM 


M-Mf  [-]■  MM 
[•!■»]  •  MT 
KKKKMMM 
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3*  POTENTIAL  ENERGY 


Tbs  total  strain  energy  of  a  multilayered  cylindrical  shell  of  radius  R  and  length 
L  can  be  expressed  as  the  cum  of  the  following  two  expressions  (Reference  17) 


n, 


,L  ,2t R 

—  f  f  (  a  NZ  +  O  N  *  + 

2  J  J  1  I!  s  22  y 


0  O 


c  N*  +  2  o  M  N 
(3$  xy  12  s  y 


♦  2o|6NxMsy  +  2a2«^y  Wsy  }  d* 


(7) 


L  2*rR 

I  f  f  .  2  ,  *  ,  £  ,  .  »  .2  .  -  , 

--  j  J  ^  d„fc„  +  2e10Ss„  fc. 


'll  “a 


28  y 


88  *y 


*12"*  y 


O  0 


♦  2d,>»  *  2d»V»»S><hd» 


(S) 


The  energy  associated  with  the  work  done  by  :he  unifr.rmly  applied  axial  com¬ 
pressive  end  load  is  given  by 

2irR  L 


*  -  /  K1  ,  d»/  u.= 

0  *-L  o 


da 


i9! 


The  total  potential  energy  is  then  given  by 
U  =  ii,  +  n2  +il3 

4.  EQUILIBRIUM  AND  COMPATIBILITY  EQUATIONS 


(IO> 


The  equilibrium  equations  resulting  .from  the  first  variation  of  the  total 
potential  energy  with  respect  to  the  displacements  u  ,  V  and  w  are  given  by 


N 

N 

N, 


♦  N  =0 

ay^y 


(ID 


+  W 

*y»*  y,y 


( w 


o 


,xx 


)  *  2N_  (M 


,xy  ,ay 


,  +  Vw.«-' 


<-  +  2M  -4-  M  =0 

a, xx  xy,xy  y,yy 


(12) 
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The  thr<5e  unknowns  in  Equation  11  can  be  replaced  by  a  single  unknown  F  (x,y), 
the  Airy  stress  function,  which  is  defined  as 


=  F 


yy 


*  y  =  F,  x  x 

I  =  -  F 
xy  ,xy 


(13) 


Utilizing  Equations  6  and  13  in  Equation  12  the  equilibrium  equation  in  radial 
direction  can  be  written  as 

di2F,xxxx  +  <  2d62-  <W  F,xxxy  +  1  dn  +  d22_  2dco)  F»xxyy 
+  (  2de|- d2e)  F  xyyy+  d2|Ftyyyy  +  Fxx(Wtyy+  -^- )  +  Fyy  W|XX 
-  2  F>xy  WXy  -•  d  n  (  ^fXKXX“W  jXRXK)  -  4d|{}(  W^xxxy-“  WjKXXy) 
(2d|S.+  4defi  )(  WjXRyy  -  W  xgyy)-  4d^(  W  syyy  -  V7  Ryyy  ) 

d2«*W,yyyy  ~  W,yyyy  *  =  0 

Eliminating  U  and  V  from  Equation  1  the  compatibility  condition  can  be 
written  as 


e  +•  e 

x,yy  y,xx 


-  « 


*yt*y 


=  (w 


»*y 


)8  -  w 


»yy 


2 

R 


xx 


~<W  „„)*  +  W 


»*y 


w 


•  »*  »yy 


—  w  Kv 

R  ,  K  X 


(13) 


The  above  equation  can  be  expressed  in  terms  of  stress  function  by  utilizing 
Equations  6  and  13  as  follows 

°22F,sxxx  ~  2a8©F,xxxy  *  {2aig+  °es)  Fxsyy  “  2ois  F,xyyy 


+  a.  F 


uF,yyyy  '  ~di2(  W,xxxx  “  W,xxxxJ  “  (  2d68~  die)(  W,xxxy_V^xxxy  1 


-{d.«+  d22-  2dee»  W,xxyy- W  xxyy)-(2d6|-d?0K 


W. 


-  W 


,*yyy 

2  }T!7  t2 


Uyyy1  d2< ( w,yyyy~  w,yyyy 5  +  (  w,*yJ  ~?w,xy} 

-  <  w  xxw  uu  -  w„  wuu)  -(  -i  w  „  -  4- w„)  (16) 


,xx  ,yy 
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5.  RADIAL  DEFLECTION  FUNCTION,  INITIA  L  IMPERFECTION, 

AND  COMPATIBLE  STRESS  FUNCTION 

The  following  trignometric  expressions  are  used  to  represent  the  total  radial 
displacement,  W  and  initial  radial  imperfection, 


W  =  W,  co#  -j-5-  co# 


-P-  ♦  W 
*y 


g  COS  -TT5-  +  Ws  COB  ZfZ-CQ*  T11 
x  x  y 


4ttx 

+  VV  COS  —B -  + 

4  ix  O 


TJ7  717  7T  X  77* 

W  -  W.  COS  -TI —  COS  T 


.  —  277-x  .  tt7  2mt  27r  y 

+  W_  COS  -3 +  W,  CO#  — 3 —  COS  —JT1- 

2  S  4Z  Jt  Cy 


where  W0  through  W4  are  the  unknown  coefficients  in  the  radial  displacement 
function  and  W,  through  are  the  amplitudes  of  the  imperfections,  12  s  and 
are  the  half-wavelengths  in  the  axial  and  circumferential  directions. 


Equation  17,  assumed  to  represent  the  total  radial  displacement,  is  of  the 
same  form  as  that  used  for  investigating  the  postbuekling  behavior  of  geometrically 
perfect  anisotropic  shells  in  Reference  17,  and  for  isotropic  shells  in  Reference  19, 
In  the  selection  of  the  buckling  configuration  (Equation  17}  assumption  is  made 
that  the  shell  may  not  deform  into  a  torsional  hackling  mode. 

Substitution  of  the  expressions  for  the  assumed  deflected  shape  and  initial 
imperfection  into  the  compatibility  equation  and  subsequent  solution  yield  the 
following  expression  for  the  stress  function 


F  =  F  -  £1- 
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In  which  cr  is  the  applied  axial  compressive  load  and  the  expression  for  F  is 
as  given  below 


•=■  _  vx  Try  .  ttk  7ry  ,  ^  2 irx  _  Zsry 

F  =  F«  cos  -jj  cos  +  F2sm  -jj-  +  Fs  cos  +  F4cos— 


.  c  _  2yx  .  2yy  .  _  .  2ttx  2try  .  ..  3vx  7ry 
+  F8C081^  cos^87+  F^.n^-sm-^  +  F^os-^  cos-^J- 

.  _  .  3trx  iry  .  _  ttx  3try  .  _  .  w*  .  3iry 

+  F8*m-|— sm-J-  +  FBcos-^--cos-|^-+  F^s.n-^sm-^- 

t  4tx  ..  4ry  .  _  4 tx  2tt  y 

+  F„cos-^-  +  P,2cos-j^-  +  Piscos-^cos-j^- 

,  _  .  4irx  .  2iry  .  5xx  Try  .  _  .  5trx  .  ,ry 

+  F^s.n-j^-  sm-^-+  ^a008  TT “* T^T  n  ~~TT 81  n  TJ* 

„  6ttx  2iry  .  _  .  6trx  .  2 iry 

+  F  coa-j —  cos  -■■- +  F.^stn  — — tm-y- 

*-K 


,-sm  — - — sin— jr 


where  the  parameters  F,  through  Fia  which  are  the  functions  of  parameters 

vs.  through  w4  assi*W,  through  "W4,  are  as  gives  below 

-{*.  ( G.  -  -frHMK}-*,6.®. 

Rr 

F  =  —g  ~  - 

6a  -  6« 


v>,®,  ♦{*,(«. 

r  2  r  2 
G2  “  G« 


_  f  d|2  *y  1  Gi 

F.  =  A.  4 - -  f  ■■  ^_2-  1 - g-  [ - ITT— 

*  2 1  o22  4^izo22Rt rz  J  32 ft  a2t 


V  c.{-  TTFTT  }  +  b»{--^T—  } 


32 17  Q„ 


217  0„ 


{  A^®''  4RW-8  )  -2BJ-A=6»«. 


r2  r2 

G  2  ~  G4 
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A’G»G*  ♦  {  *>(«.-  )  -2bJg4 


6  - 

-®«2 

II 

ll7 

-  {4B3  «•  2B, 

+  eS4] 

K '<<*.* 

-<*> 

F8  = 

{4B,  +  2B,  + 

81!«}  <V<GS2- 

G«2) 

F9  = 

-4  83G7/(Gt2- 

^0  = 

4®3G#/(Gr2-G 

;> 

F(|  = 

A.  + 

V 

> 

1 

G3 

II 

4  1  °22 

,6^2°2^'2r  ^ 

32^  '«!, 

FI2  = 

_-n 

Ul 

II 

-ff.G./<e,8- 

F  = 

14 

B  <5  '(G  a-6  2) 

8  10  9  10 

FI8  = 

-  8  G|(  f  (  G|f2- 

G,l> 

F»e  = 

8«4G,2^6,|2-- 

G«> 

f,t  = 

-2B2G,3/(G,s2 

"  Gt4  » 

F«  = 

.  2, 

,,4> 

tl 

14. 

S.  y 

V  = 

5* 

17 

A.  =  W.  -  ",  a2  =  wg  -  W2 

A5  =  W3  ~  Ws  -  W4 
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9t 

*"3 


W,  *3 


w  w 

I  3 


B 

4 


w  w  -  w  w 

I  4  14 


G.  =  d2,,72+  (d„  +  d22-2d63) 

G2  8  *  °|.’?2+  <2a|a  +  06e) 

Gs  ^  2dfl2  -  d,a )  /i.  +  {  2  d4J  -  d£6)  77 

°4  8  2/Ia2«  +  2^0,6 

g8s  ®,°22^8+  °„^2  ♦  9(2o,2+  aee) 

Ge  s  54/r°26+  6^°,« 

G7  =  *  9lanV*  +9(20,2  *  °66i 

Ga  s  6?a26  +  54  ’?°|« 

6,=  Ma^/T*  ♦  20,,^*+  8(2o,2  +  a^) 

°»*  *2Pen*  8^a.o 

G„=  625  Ogg/T  2  ♦  o,,^3  -  25  (  2g|?  ♦  <s@#> 

G,2=  250/t  ags  ♦  «0^o{e 

G  =  BlflJI1  +  «»*♦»(  2a  -«■  a  ) 

a  22^  15  '  !2  ©8 

G  =  54  a  IT  +  6  a  Tj 

14  28“  18  » 
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6.  UNIT  END-SHORTENING  AND  UNIFORM  RADIAL  DISPLACEMENT 


The  unit  end-shortening  caused  by  the  uniformly  applied  compressive  end  load 
is  given  by 

L 


e  *  -f-  /  u.> 

0 


dx 


(22) 


then,  through  utilization  of  Equations  1,  17,  18,  and  19,  and  integrating  Equation  22 
end-shortening  can  be  written  as 


£  =  au€r  +[(W12-W*)  ♦  8(W22-W22! 

♦  4<  W  2-  W  *)  «•  32 (  W  2—  W2'|  -- 

9  9  4  4  J  at 


(23) 


The  continuity  condition  for  the  circumferential  displacement  is  given  by 
2irR  2ttR 

'  V,y  d>'  =  *  °22F,,X  -  °»F»,  +  dl2W,ltK+  V,„ 


+  2de2'!',xy-Tlw.y,2tT(w,/+f--T]!lv  =0  ,2',) 


Substitution  of  the  assumed  radial  deflection  function,  initial  imperfection  and 
stress  function  into  the  continuity  condition  yields  the  following  expression 
for  W0 


%  =  ■F[(w,J-vy,ai  +  <kw|-w5!iJ 


2  —2x1  ^  R 

8iJ 


(25) 
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7.  TOTAL  POTENTIAL  ENERGY 


Utilizing  Equations  17  through  25,  the  following  expression  for  the  total 
potential  energy  is  obtained 

=  TX  +  IBAjj+AT.A.B,/ 

+  64  ViV’*  4B>*  +  64  Bl  /**}  +  -rf-pr  {  2  V? 

‘l  *2 

♦  32  V !  +  4  *,  B,M  +  64  A3  B2  /  }  ♦  — -X,  {  (  A  B, 

8  '"'a 

+  zb,  +  ea4 12  ♦  6«b|}  /  *  -rht  '«b. 3V 

•t  "*  *  n 


where 


Te  „  .  „  2  4  Tn 

- F  64B4^  +  ~Zz~- 2 


(9  T)0 


r.  -  m2 


4  B2yu.4  +  2A2oj2  +  2  Aftu2 


♦  IZ  c,2  +  T  C3  -  T  C,  V  -  2C3a4“  +  <  *  16 

4-{32A2V2+5IZAaK24sLc3+  -jpC2  +  SB2  4  C,B„ 

+  Af<^>  +  32  a|<^>  +  i6A^<£>  +  5 1 2  A  2<£  }  ~k\  A2+  4  A2 )  I6uyi2 

+  «-(AzC(  +  l6A4Csi2/i2-*-  />V^T  (A2  +  32  A2 )  2  \x 
-  {C,  +  8C2+  4C3+  32C4|  2cr/u.2tu|  +  constants  (26) 


P.  = 

i 


Wj 


(°22d22) 


Pi  = 


W, 


.  .a**1'2 

(°22  d22J 


i  =  I  through  4 


/*  = 


N  = 


lz_ 

7rR 


(^f 


u>  = 


,n  J/2  w2 
(°22d22)  N 


(27) 


er  R 


(  -t 4-  r 
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A,  =  P,  -  F> 

A2  =  P2  -  Pg 

A3  =  P3  -P3 

\  =  P4'P4 

8,  =  P.P2  -  P,P2 

B2  =  P3P4  "  P3P4 

B  =  PP  -PP 

3  13  13 

B  =  P  P  -  PP 

4  14  14 

05=  P2P3-P2P3 

c,  =  P  2  -  p" 2 

2  — 2 

C2  =  P2  "  P2 

C  =  p  2  -  p"  2 

“3  a  KS 

C  -  P  2  -  F  2 

*  ^4  K4 

The  values  of  T(  through  T|3 

appearing  In  the  energy  expression  are  as  follows 

T,  = 

4  2 

/l  +  ClfJL  +  1 

T0  =  6/i  (/i2/3  +  9y ) 

ll 

2/1  (/i2yS  ♦  y  ) 

T9  =  625 /i4  +  25  a/i2  +  1 

T3  » 

K /l*  +  1^/12+  X  - 

/i2a>  T|0  =  lO/i  ( 25/t2/3  +  y) 

= 

/i(  v/i2  +  £) 

Tn  =  I6/x4+  4  a/i2  +  1 

T*  = 

Si  /i4  +  9  a/i2  +  I 

r|2  =  4/i  (  4/i2£  +  y  ) 

II 

»-* 

6xi  (  9/i2/S  +  y  ) 

TI3  =  +  V'M2  +  X  ~  -y-  ^ 

It 

/i4  +  9  e/i2  +  i 

where  the  stiffness  parameters  a, (3  ,  etc.  are  defined  by  the  following  relations 

a 

2oi2+  oss 

-/  Q  Q 

V  «  1  *  AA 

1  t  C.C. 

°  26 

S/4  3. '4 

a..  a__ 

:*  -  t  i. 
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16 


n  1/4  „  3,4 
°22  Qll 


c  = 


K  - 


A  = 


12  /j 
22  V 


Oil 

d 


22 


'2! 


■ii 


+  d22  -  2d 


69 


J  °ll  *tz 


*12 


+  2d 


66 


(28) 


*  = 


'll  q  it 


‘22  u22 


The  end-shortening  6  (Equation  23)  may  be  redefined  as  a  new  parameter 
6  =  *  +  -kK-'T>  +  8IS>!ii-Ps2)  MIP'-P,*! 

+  32(P/-P42)]-^  (29) 


where 


6  = 


SR 


2  Ald22 


(30) 


14 


AFFDL-TR-68-13S 


8.  NONLINEAR  ALGEBRAIC  EQUATIONS 


The  condition  that  for  a  given  c r  (i.  e. ,  applied  axial  compressive  load)  the 
total  potential  energy  has  a  stationary  value  for  small  variation  of  P  ,  P  ,  P3 ,  P 
is  equivalent  to 


dv  _  ^  _  t?V  _  dy_ 
dPt  =  dP2  =  dP3  "  dP4 


=  0 


(31) 


This  yields  the  following  four  nonlinear  algebraic  equations 

jr^r{z{T?*T*'*,+  4VB,  +  a,ps)M2+  bb.p^4} 


■l  '2 


rjfr{«V.  +  V<b,  +  a,p3>}  ♦  {  ziab^sb, 


JjT« 

T2-  . 

I  2 


5  '6 


♦  8B4)(4Ps  +  2Pa  ♦  8P4  )^4}  +  t/T  .  *2  B3PjM‘ 

*7  “  T$ 


♦  7;  T-„7  123  8^^^'*  4-  S^._F»aatu  +  2A,  +  { 


T  -  T 
'9  1 10 


C,P, 

-V  +  2F?Ba 


+  2A,<£  |  fi*  +  4  kA2P(  yti2  +  4/3  V^T  A^2  -  4P  <r/i2  w  -  0 


T, 


C|  o> 
2 


+  2«C{fi 


{4T3A,P|M2+  8B,P,/}  +  yV 

'l  ”  '2 

-  {  8 (J( 2 B3+  8,  +  4B4 )  flA  }  +  8B.P,/  +  4A3«j2 

Tll  TI2 

{64  A2x2+  ISBsP3  +  C(P3  +  64  A2^>}^4  -32/cA2o)^2 
2  -  32  P2  crp  cj  =  0 
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■prjt  {32A,(T;J+r,aM-S-n;j/(Ba+  A,P,)  +  128  8,?,/} 
■’■TT^T  {64t,jA3+  64M2(B2+  A,?,)}  +  J\y  {  I6(2B5+  B 

't~  *9  '5  ~  *6 


+  «B,'p,  +  '“Vi}/*  zrbi32^^  +  t«T|' 8  8B.1M4 


T  -  T 
'r  1  s 


T  -  T 
'll  12 


+  64xA4P3M2+  2CsP3- 4A4P.o)  +  32  A3  +  {  2  Cs  P3  +  16  B5  P2  +  C,  P2 


+  32  j  ju4  +  128  A^2  ~  ISPc/t' 


w  =  0 


^7?  {«4T,3A  3P3M2  +  128  S2  (j  p.4  }  -3-  ?T<  ,  ««A,P,M' 

■  »*“•«  1*2 


-2 


+  -ayia  {  32(2B3+  B,  +  4B4)F>  +  I26B2P5}/ 
*8  'e 


+  — ir~ — g-  *28  8^  P  fi*  +  4fi4u2  -  2 CjOJ  +  {  #c2+  1024 

To  Ta 


s  *10 


-  128  k&4  ujfi2  +  32  #cC3/tt2  -  !28  o* ^t20i  F^(  =  0 
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The  iir perfection  sensitivity  of  a  cylindrical  shell  can  be  measured  by  a 
quantity  p  *  which  is  defined  as  the  ratio  of  the  buckling  load  of  an  initially  im¬ 
perfect  shell  to  the  classical  buckling  load.  If  <7m(JX  is  consioered  as  the  buckling 
load  of  the  imperfect  shell  and  cr  its  classical  buckling  load  then 


/>*  = 


max 


Tcl 


(32) 


For  small  imperfections  the  cylindrical  shell  passes  through  a  maximum  load 
crmQX  state  before  it  snaps  through.  In  case  of  large  imperfections  however,  the 
cylinder  does  not  experience  snap-through  and  hense  crmQxis  not  well  defined 
(see  Figure  2),  For  a  given  shell  crmQ%  depends  upon  imperfection  parameters 
P,  through  P4 ,  p.  ,  and  w  . 


The  values  of  deflection  parameters  P,  through  P4  can  be  obtained  by  solving 
Equation  32.  The  solution  of  this  set  of  equations  depends  upon  the  parameters  a 
through  (stiffness),  P(  through  P4  (initial  imperfection  amplitudes),  p.  (aspect 
ratio),  at  (circumferential  wave  number)  and  <y  (applied  axial  load).  End¬ 
shortening  £  corresponding  to  the  applied  load  c r  can  be  calculated  from 
Equation  29.  The  stiffness  parameters  a  through  <f>  can  be  calculated  knowing 
the  thickness  of  the  shell,  number  of  layers  and  the  elastic  properties  of  material 
in  each  layer  (Reference  17).  P,  through  P4  depend  upon  the  amplitudes  of  as¬ 
sumed  deviations  from  the  geometrically  perfect  cylindrical  shell.  The  parameters 
p.  and  w  depend  upon  the  wavelengths  of  assumed  initial  georr  -trie  imperfection. 
For  isotropic  cylindrical  shells  Koiter  (Reference  20),  Hutchinson  (Reference  21), 
and  Madsen  and  Hoff  (Reference  8)  predicted  that  the  imperfections  with  wave¬ 
lengths  equal  to  those  for  small-deflection  buckling  mode  are  most  critical.  The 
parameters  p  and  oj  for  anisotropic  cylindrical  shells  corresponding  to  buck¬ 
ling  mode  for  small-deflection  theory  can  be  calculated  from  equations  given  in 
Reference  17. 


9.  NUMERICAL  ANALYSIS 

The  set  of  nonlinear  algebraic  equation0  (Equation  32)  can  be  solved  by  Newton- 
Raphson  iterative  technique  (Reference  22).  A  FORTRAN  IV  computer  program 
was  written  for  IBM  7094  to  solve  the  equations  and  to  evaluate  £•  for  given  cr  , 
or  c  for  given  £ . 
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SECTION  in 

RESULTS  AND  CONCLUSIONS 


From  the  theoretical  analysis  developed  in  the  previous  sections,  it  is  seen 
that  the  behavior  of  the  cylindrical  shell  depends  upon  ten  stiffness  parameters 
a  ,  /3  ,  etc.  Because  of  this  large  number  of  parameters,  it  Is  not  possible 
to  present  meaniful  general  results  in  terms  of  stiffness  parameters.  Accordingly, 
a  three-layer  cylindrical  shell  of  radius  6. 0  inches  and  thickness  0. 036  inch  is 
selected  to  illustrate  the  application  of  the  theory  and  to  arrive  at  conclusive  results. 
The  thickness  of  each  layer  is  0.012  inch  and  all  layers  consist  of  either  giass- 
epoccy  or  boron-epoxy  composites.  The  elastic  constants  of  the  composites  are 
as  given  bslow. 


glass  -  epoxy 

E„  =  7.5  x  IOS  psi 
£3,,,=  3.5  x  10*  psi 

,l2=  0.25 

6  =  1.25  x  I06  psi 

"a. 

Two  sets  of  fiber  orientations  are  considered.  The  fiber  orientations  in 
outer,  central,  and  inner  layers,  respectively,  are  I )9?-9°,0°,  2)  9°,  "6°,  90° 
where  6  assumes  values  ranging  from  0°  to  90°  with  10-degree  intervals.  The 
two  sets  will  be  referred  by  notation  1  ■  9  and  2*6,  and  the  numerical  value  of 
9  will  indicate  corresponding  fiber  orientation.  For  example,  1. 20  corresponds 
to  fiber  orientation  <20°,  -20°,  0°).  When  $  is  zero,  fibers  are  axially  oriented, 
and  when  it  is  S0°  the  fibers  are  circumferentially  oriented  (Figure  4).  The  stiff¬ 
ness  parameters  3  ,  0  p  vtC*  §  for  these  sets  of  fiber  orientation  may  be 
found  in  Reference  17.  The  critical  classical  buckling  load  parameter  crC£  , 
buckling  load  er=  ,  aspect  ratio  parameter  p.  and  the  circumferential  wave 
number  for  small-deflection  theory  are  given  in  Tables  1  and  2  for  glass-epoxy 
and  boron-epoxy  composite  cylindrical  shells.  These  values  are  calculated  with 
the  aid  of  the  computer  program  given  in  Reference  17. 


boron -epoxy 


22 

*.2  8 


=  vir,  E/  E 
12  22  11 


40.0  x  10  psi 

4.5  x  10°  psi 
0.25 

1.5  x  !0S  psi 
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To  study  the  influence  of  initial  geometric  imperfections,  the  relative  values 
of  the  initial  imperfection  amplitudes  and  parameters  are  set  to 

Wj  =  +4^  =  -4¥3  =  16  W4 

and  _____ 

P,  =  *4P2  =  -4p3  =  ,6P4 


These  are  the  same  as  those  given  by  Yoshimura  buckling  pattern  (Reference  8). 
The  imperfection  ratio  W  *  is  defined  here  as  the  ratio  of  the  sum  of  initial  im¬ 
perfection  amplitudes  to  the  thickness  of  the  shell  t  .  The  following  relation  is 
obtained  from  Equations  27  find  34. 


W  *t 


For  a  specified  value  of  W  *  the  initial  imperfection  amplitude  parameters  F>  . 
through  P4  can  be  calculated  from  Equations  34  and  35. 


For  isotropic  cylindrical  shells  all  stiffness  parameters  given  in  Equation  27 
are  zero  except  that  a/2  =  P  =  <£  =  »,  Figure  3  shows  the  curves  of  the  end- 
shortening  parameter  £•  against  the  axial  load  parameter  cr  for  an  isotropic 
shell  for  different  values  of  imperfection  ratio  W  * .  When  W*  is  zero  crft)aa  is 
equal  to  <rc ^  which  is  equal  to  unity  for  the  isotropic  shell.  In  Reference  8, 
only  two  terms,  corresponding  to  W,  and  "w2  are  used  in  the  assumed  W*  function 
to  study  the  behavior  of  an  isotropic  shell.  It  is  observed  that  for  an  isotropic 
shell  the  values  of  crmaK  evaluated  from  the  theory  developed  for  anisotropic 
shells  in  this  report  are  comparable  to  those  given  in  Reference  8. 


Figures  4  and  5  show  the  load-shortening  curves  for  glass-epoxy  and  boron- 
epoxy  composite  cylindrical  shells  for  fiber  orientation  (0°,  0°,  0°),  In  Tables 
3  to  18,  the  values  of  <7maa ,  'cFmax  P*  and  &mox  are  given  for  W  *  =  0. 01, 
0.02,  0.04,  0.06,  0.08,  0.1,  0.2,  and  0. 3  for  glass -epoxy  and  boren-epoxy 
composite  cylindrical  shells.  The  same  results  are  also  presented  in  Figures  6 
through  13.  The  discontinuity  in  the  curves  simply  indicate  that  for  that  range 
of  orientation  there  is  no  snap-through.  Figure  6  gives  the  variation  of  <rmax 
with  the  fiber  orientation  (  Q°  ~6°  0°)  for  a  three-layer,  glass-epoxy  com¬ 
posite.  Figure  7  shows  the  variation  of  imperfection  sensitivity  p*  for  the  same 
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set.  From  these  figures,  it  f  3  evident  that  increase  or  decrease,  in  classical 
budding  load  with  change  in  fiber  orientation  is  generally  accompanied  by  a  de¬ 
crease  or  increase  in  imperfection  sensitivity.  The  same  behavior  can  be  ob¬ 
served  in  the  case  of  a  glass-epoxy  composite  for  fiber  orientation  (  6 ,  -6,  90°). 
From  Figures  8,  9,  12,  and  13  similar  behavior  is  observed  for  boron-epoxy 
composite  shells.  In  Figures  7,  9,  11,  and  13  the  values  of  a*  for  isotropic 
shells  are  given  for  the  purpose  of  comparison.  It  can  be  seen  that  the  composite 
shells  are  less  inperfection  sensitive  than  isotropic  shells  and  boron-epoxy 
composite  shells  are  .less  imperfection  sensitive  than  glass-epoxy  composite 
shells. 
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Figure  i.  Geometry 


Figure  4.  Initial  Imperfection  Solution  of  Glass -Epoxy  Composite  Cylinder  for  Fiber  Orientation  (0 
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FABLE  1. 

CLASSICAL  BUCKLING  LOADS  FOR  GLAS5-6PQXY 
COMPOSITE  CYLINDRICAL  SHELL 


CASE 

1/ 

r~ 

N 

I.  C 

1 . 58Q8E-0 1 

4.9158E 

02 

1.C330E  00 

11 

1,10 

7 .81546-0 1 

4.99376 

02 

9.97406-31 

12 

1  »  2v 

8.4265E-01 

5,15996 

32 

9.32536-01 

12 

i  o  30 

9  «  22556-01 

5.3507F 

32 

9,30276-01 

12 

1.43 

9.785  36-01 

5.4473E 

02 

8.43746-01 

12 

1.53 

9  *  ?3c  46-01 

5.33746 

02 

7.9899E-31 

12 

1.60 

9. 16456-01 

5.22226 

02 

8. 03126-01 

12 

*“»• 

0 

O 

8.42306-01 

5.02956 

02 

8.32606-C1 

12 

1.83 

7.92276-01 

4.8753E 

0  2 

7.86636-01 

11 

1.90 

7.37486-01 

4.80556 

02 

8.38546-01 

11 

2.  0 

6. 87896-31 

4.80546 

02 

1.23686  00 

11 

2.10 

f • 13806—31 

4.87506 

02 

1.2712E  03 

11 

2.20 

7.71716-31 

5,02956 

02 

1.201 16  00 

12 

2.30 

8 .57896-01 

5.22226 

02 

1.24526  00 

12 

2.43 

9.361 7E-0 l 

5.3874E 

02 

1,25166  00 

12 

2,50 

9.60706-31 

5.44736 

02 

1.18526  00 

12 

2  •  6u 

9.1760E-01 

5.35C7E 

02 

1.07496  00 

12 

2.70 

S.439QE-01 

5.  16C3E 

02 

1.31776  00 

12 

2.80 

?. 82726-01 

4.99066 

02 

1 . 0026t:  00 

12 

2.9C 

7.58386-01 

4.9158E 

02 

1.00006  00 

11 
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TABLE  2. 

CLASSICAL  BUCKLING  LOADS  FOR  BORON-EPCXY 
COMPOSITE  CYLINDRICAL  SHELL 


CASE 

aC& 

**1 

W 

1.  0 

4.9392E-01 

8.2929E 

02 

l.OOOOt  00 

10 

i .  ic 

5.6225E-01 

9.2130E 

02 

1.0012E  00 

11 

1.20 

7.3151E-01 

1.1134E 

03 

9.7972L-01 

12 

1.30 

9.4227E-01 

1.2958F 

03 

8 . 7935E-01 

13 

1.40 

1.0168E  00 

1.3212E 

03 

7.4584E-01 

13 

1.50 

9.1515E-01 

1.2163E 

03 

7.2442E-01 

13 

1.60 

7.6007E-01 

1.0883E 

© 

03 

7.9051L-01 

13 

1.70 

6.0374E-01 

9.6170E 

02 

8.9527E— 01 

13 

1  .BO 

4.5867E— 01 

8.4377E 

02 

9.4702E-01 

12 

1.90 

3.7875E-01 

7.7412E 

02 

9.6256L-01 

11 

2.  0 

3.6541E-01 

7.7412E 

02 

1.0389E  00 

11 

2.10 

4.1 170E-0 i 

8.4378E 

02 

1.0559E  00 

12 

2.20 

5. 1807E-01 

9.6170E 

02 

1.1170L  00 

13 

2.3u 

6.8294E-0I 

1.C883E 

03 

1.2650E  OC 

13 

2.40 

8.8459E-01 

1.2164E 

03 

1.J8C4E  00 

13 

2.50 

1.0151E  00 

1.3212E 

03 

I.3407E  00 

13 

2.60 

9.5486E-01 

1.2958E 

03 

1.1372L  00 

13 

2,70 

7.5876E-01 

1.1133E 

03 

1.0207E  00 

12 

2.80 

5. 7806E-0 1 

9.2131E 

02 

9.9877E-01 

11 

2.90 

4.9392E-01 

B.2929E 

02 

l.OOOOE  00 

10 
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TABLE  3. 

BUCKLING  LOADS  FOR  IMPERFECT  GLASS-EPOXY 
COMPOSITE  CYLINDRICAL  SHELL  FOR 
_*■ 

W  =  O.Ol 


CASE 

O'mas 

P* 

£max 

i.  0 

6.8687E-01 

4.4540E 

02 

9.0607E-01 

6.9242E— 01 

i.10 

7.0550E-01 

4.5051E 

02 

9.0270E-01 

7 . 1 192E-01 

1.20 

7.5472E-0I 

4.621SE 

02 

8.9565E-01 

7.5992E-01 

1.30 

8.1979E-01 

4.7547E 

02 

8 . 886  IE-01 

8.2392E— Cl 

1.40 

8.70486-01 

4.8438E 

02 

8 . 8922E-01 

8. 7392E-01 

1.50 

8=721 7E-01 

4.8249E 

02 

8.9560E-01 

8.7544E-01 

1.60 

8.2606E-01 

4.7071E 

02 

9.0137E-01 

8.2994E-01 

1.70 

7.6338E-C1 

4.5583E 

02 

9.0630E-01 

7.6814E-01 

1.80 

7 .16796-01 

4.4669F 

02 

9.1629E-01 

7.2098E-01 

1.90 

6.9466E-01 

4.4071E 

02 

9.1709E-01 

6.9942E-01 

2.  0 

6.2799E— 01 

4.3870E 

02 

9 . 1292E-01 

6. 3844E-G! 

2.10 

6.4861E-01 

4.4485E 

02 

9.1251 E-0 1 

6.5969E-01 

2.20 

6.9812E-01 

4 . 5499E 

02 

9 . 0464E-0 1 

7.08i2E-01 

2.30 

7 . 7536E-0 1 

4.7198E 

02 

9.0380E-0I 

7.8594E-01 

2.40 

8.4515E-01 

4.8636E 

02 

9.0277E-01 

8.5562E-01 

2.50 

8.6120E-01 

4.88 3  IE 

02 

8.96**E-0I 

8  =  S942E~01 

2.60 

8 . 1 782E-0 l 

A.7689F 

02 

8.9I26E  01 

8.241 TF-Oi 

2.70 

7.5655E-01 

'*  -62  59F 

02 

6.9649F-01 

?=6?h?E-01 

2.80 

7.3679E-01 

4.5C65E 

02 

9.0299E  di 

7. 1 342E- 01 

2.90 

6,8 01 

'•  '*'i  l3f 

1? 

9  s  0596F-0 1 

£*9J9/fc  cl 
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TABLE  4. 

BUCKLING  LOADS  FOR  IMPERFECT  GLASS-EPOXY 


CASE 

COMPOSITE  CYLINDRICAL  ShELL  FOR 
—a 

W  =  0.02 

^mox  ^ ,,0^  p » 

1.  0 

6.49006-01 

4.2085E  02 

8.56116-01 

6.59346— 01 

L.10 

6.6603E-01 

4.2531E  02 

8.5220E-01 

6.78126-01 

1.20 

7  •  1020E-0 1 

4.3489E  02 

S.4282E-01 

2. 09206-02 

1.30 

7.6875E-01 

4.4587F  02 

8 . 3329E-QI 

7.77426-01 

1.40 

8.1533E-01 

4 .53696  02 

8.3288E-01 

8.21926-01 

1.50 

8.1846E-01 

4.52786  02 

8.4045E-01 

8. 2492 E— 01 

1.60 

?. 7707E-01 

4.4280E  02 

8.4791E-01 

7.8392E-01 

1.70 

7.2Q04E-01 

4.29956  02 

8.5485E— 01 

7.28426-01 

I 

6.79366-01 

4.2337E  02 

8. 6845E-01 

6.8676E-01 

1  .^0 

Cl 

4.1796E  02 

8.69756-01 

6.66426-01 

2  •  0 

5.9681E-0I 

4.1691E  02 

8.6763E-01 

6.1687E-01 

2.13 

6. I65IE-0I 

4.2263E  02 

8.67  356' Oi 

6. 37816-01 

"3  t  a 

6.61 346-01 

4.3IG2E  02 

8.56986-01 

6.81566-01 

2.30 

7 . 3459E-0 I 

4.47166  02 

8 .56256-0 1 

7 . 5594E-01 

2.40 

B.0C35E-01 

4.6058E  02 

8.54926-01 

8.21 566-01 

2.50 

8 . 12956-01 

4.60956  0? 

8 . 462 16-0 1 

8. 3000E-Q1 

2.60 

/»6945t‘-f)f 

4.4868E  02 

5.3855E-01 

f .82196- 01 

2.70 

7. I256E-01 

4.35696  32 

8.44376-01 

7. 243ft- 01 

2.63 

6.6  ft  Jl  0  J 

4»2543f  02 

8 »  52 6t> i  —0  i 

&.7«n9£-0I 

2.90 

6  *49066-01 

4.  ?.  cast  o? 

8 . 56  L 1 L  01 

fc . 5 7346- 01 
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TABLE  5. 


BUCKLING  LOADS  FOR  IMPERFECT  GLASS-EPOXY 
COMPOSITE  CYLINDRICAL  SHELL  FOR 

H*  =  0.04 


CASE 

G’mas 

^max 

P 

1  •  0 

5.9523E-01 

3.8598E 

02 

7 . 85166-01 

6. 1562E-01 

1.10 

6.1091E-01 

3.90UE 

02 

7.8I67E-01 

6.3331E-01 

1.20 

6.4857E-01 

3.9715E 

02 

7.6968E-G1 

6.7000E-01 

1.30 

6.9807E-01 

4.0487E 

02 

7.56676-01 

7. 1341E-01 

1.40 

7.3874E-01 

4.1108E 

02 

7.5464E-01 

7.5192E-01 

1.50 

7.4257E-01 

4.1080E 

02 

7.6252E-01 

7.5492E-01 

1.60 

7.0698E-01 

4.0286E 

02 

7. 7143E-01 

7.2C42E-01 

1.70 

6. 5750F-01 

3.9260E 

02 

7.8060E-01 

6. 7384E-01 

1.80 

6.2365E-01 

3.8865E 

02 

7.9723E-01 

6.3692E-01 

1.90 

6.0549E-01 

3.8413E 

02 

7. 99356-01 

6 . 1984E-01 

2.  0 

5.5331E-01 

3.8653E 

02 

8.0436E-01 

5.9469E-01 

2.10 

5.7194E— 01 

3.9226E 

02 

8.0464F-01 

6.1656E-01 

2.20 

6. 1125E-01 

3.9837E 

02 

7 . 920  7E-01 

6. 54696-01 

2.30 

6.7934E-01 

4.1353E 

02 

7.9187E-01 

7.26256-01 

2.40 

7.3961E-C1 

4.2563E 

02 

7.9004E-0I 

7.8531E-01 

2.50 

7.4753E-C  l 

4.2386E 

02 

7.781 1E-0 l 

7 . 8375E-0 1 

2.60 

7.0363E-01 

4.1030E 

02 

7 . 6682F-0 1 

7 . 3000E-01 

2. 70 

6.5I87E-CI 

3.9858E 

C2 

7 . 7245E-0 1 

6.75946-01 

2.80 

6.1238E-01 

3.9C45E 

02 

7.8237E- 01 

6 . 3750E-01 

2.90 

5.9523E-Q1 

3.8598E 

02 

7. 651 8E-01 

6.1562E-01 
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TABLE  6. 

BUCKLING  LOADS  FOR  IMPERFECT  GLASS-tPUXY 
COMPOSITE  CYLINDRICAL  SHELL  FOR 

0.06 


CASE 

°fnox 

°max 

P  * 

^mait 

1.  0 

5.5524E-01 

3.6005E 

02 

7.32436-01 

5.8656E-01 

1.10 

5.7060E-01 

3.6437E 

02 

7 . 3010E-01 

6.0906E-01 

1.20 

6.0370E-01 

3.69676 

02 

7 . 1643E-0 1 

6.3750E-01 

1.30 

6.4738E-01 

3.7547E 

02 

7.01736-01 

6.7437E-01 

1.40 

6.8234E -01 

3. 79c9E 

02 

6.97036-01 

7.05 16E-01 

1.50 

6.8640E— 01 

3.7972E 

02 

7.0484E-01 

7.3534E-01 

1.60 

6.5463E-01 

3.73C36 

02 

7 . 143 IE-Q 1 

6. 7500E-01 

1.70 

6.1055E-C1 

3.6457E 

02 

7.24866-01 

6.35006-01 

1.80 

5.8062E-01 

3.61836 

02 

7 . 4222E-01 

6.0000E-01 

1.90 

5 .64296-0 1 

3.5799E 

02 

7.4496E-01 

5.8531E-0L 

2.  0 

5.2161E-01 

3.6438E 

C2 

7.58286-01 

5.8812E-01 

2.10 

5.3961E-01 

3.70C9E 

02 

7.59166-01 

6. 12 19E-01 

2.20 

5.7570E-01 

3.7520E 

02 

7.4601E-01 

6.4875E-01 

2.30 

6.4032E-01 

3.8978E 

02 

7.46396-01 

7.2000E-01 

2.40 

6.9659E-01 

4.0087E 

C2 

7.4408E-01 

7 . 7344E-01 

2.50 

7.0083E-01 

3.9738E 

02 

7.2950E-01 

7.5969E-01 

? .  60 

6.56306-01 

3.8270C 

02 

7.15246-01 

6.9844E-0L 

2.7  0 

6.C783E-01 

3.7166E 

02 

7. 20266-01 

6.45946-01 

2.80 

5.721 76-01 

3.6481E 

02 

7.310:E-0i 

6.1156E-01 

2.90 

5. 55246-01 

3.6J-5E 

02 

7. 32436-01 

5.86566-01 
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TABLE  7. 

BUCKLING  LOAOS  FCR  IMPERFECT  GLASS-EPOXY 
COMPOSITE  CYLINDRICAL  SHELL  FOR 

¥*=  0*08 


CASE 

a 

max 

a 

max 

P* 

^max 

l.  0 

5.2289E-01 

3.3907E 

02 

6.8976E-01 

5.6562E-01 

l.io 

5. 3849E-01 

3.4386E 

02 

6.8901E-01 

5.9281E-01 

1.20 

5.6800E— 01 

3.4781E 

02 

6.740&E-Q1 

6. 1563E-01 

1.30 

6.C669E-01 

3.5187E 

02 

6.576ZE-01 

6.4406E— 01 

1.40 

6.3786E-01 

3.5494E 

02 

&.5159E-0L 

6.6625E-01 

\  .  50 

6.4115E-01 

3.5469E 

02 

6. 5837F-01 

6.6687E— 01 

1.60 

6. 1227E-01 

3.4889E 

02 

6.6809E— 01 

6. 40 3 IE- 01 

1.70 

5.7246E-01 

3.4182E 

02 

6.7964E-0! 

6.0562E-01 

1.80 

5.4496E-01 

3.3961E 

02 

6.9664E-01 

5. 7063E-01 

1.90 

5.3016E-01 

3.3634E 

02 

6 . 9990F-0 1 

5.5812E-01 

2.  0 

4.9649E-01 

3.4683E 

02 

7.2 1 76E-Q1 

5.93448- 01 

2.10 

5. 14138-01 

3.5261E 

02 

1 .2331E-01 

6. 2094E-01 

2.20 

5.4842E-01 

3.5/42F 

0? 

7. 10668-01 

6.62 19E-01 

2.30 

6.1059E-01 

3.7168E 

02 

7 • l 1 73E-0 1 

7  s  3687E-01 

2.40 

6 . 6360E-0 l 

3.-8188E 

02 

7.0805E-01 

/ .8312F-01 

2.50 

6.6443F-01 

3.7674F 

02 

6. 91618-01 

7.51 56E-01 

2.60 

6. 18998-01 

3 .6094 E 

02 

6.7457E-01 

6.7937E-01 

2.70 

5 . 7293F-0 1 

3.5030E 

02 

6.788IC  Oi 

6.26878-01 

2. 80 

5.401 6E-01 

3.4440E 

02 

6. 9C1 IE-01 

5.9594E-0i 

2 .90 

5. 22898-01 

3.390 IL 

0? 

6.897b8-01 

5.64628-01 
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TABLE  8. 


BUCKLING  LOADS  FOR  IMPERFECT  GLASS-tPOXY 
COMPOSITE  CYLINDRICAL  SHELL  FOR 

W*=  O.IO 


CASE 

°max 

^max 

P* 

^mox 

1.  0 

4.9563E-01 

3.2139E 

02 

6 . 5380E-01 

5. 5094E-01 

1.1C 

5.1189E-01 

3.2688E 

02 

6.5498E-01 

5.8469E-01 

1.20 

5.3838E-01 

3.2967E 

02 

6.3891E-01 

6.0156E-0A 

1.30 

5.7282E-01 

3.3223E 

02 

6.2091E-01 

6.2187E-01 

1.40 

6.0037E-01 

3.34006 

02 

6.1329E-01 

6.37196-01 

1.50 

6.0317E-01 

3 .33666 

02 

6. 1937E-01 

6.3625E-01 

1.60 

5.7660E— 01 

3.2856E 

02 

6.2917E-01 

6. 1250E-01 

1.70 

5.4036E-01 

3.22666 

02 

6. 41536-01 

5.8312E-01 

1.60 

5.14366-01 

3.2054E 

02 

6.5752E-01 

5.4656E-01 

1.90 

5.0089E-G 1 

3.1777E 

02 

6.6126E-01 

5.35946-01 

2  •  0 

4. 75836-01 

3.3240E 

02 

6. 91 72E-01 

6.1031E-01 

2.10 

4.9333E-01 

3 .38356 

02 

6.9405E-01 

6.4375E-01 

2.20 

5.2730E-01 

3.4366F 

C2 

6.8329E-01 

7 • 1875E-01 

2.30 

2.40 

6.3811E-01 

3.6721E 

02 

6.8162E-01 

8.4156E-01 

2.50 

6. 3505E-0 1 

3.6008E 

02 

6.6103E-01 

7 . 5969E-01 

2.60 

5.8832E-01 

3.4306E 

02 

6.41 15E-01 

6. 7031E-01 

2.70 

5. 44046-01 

3.3265E 

02 

6.4467E-01 

6.1656E-01 

2.80 

5.1366E-01 

3.2751E 

02 

6.56256-01 

5.8844E-01 

2.90 

4. 95636-01 

3.2139E 

02 

6. 5380C-01 

5. 5094E-01 

**  NC  SNAP-THRCUGH 
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TABLE  9. 

BUCKLING  LOADS  FOR  IMPERFECT  GLASS-bPOXY 
COMPOSITE  CYLINDRICAL  SHELL  FOR 

W*=  0.20 


CASE 

<T 

max 

< T 

max 

P 

^max 

1.  0 

4.C216E-C1 

2.6078E  C2 

5.3C50E-01 

5.4406E-01 

I.  10 

** 

1.20 

$$ 

1.30 

4.61 52E-0 1 

2.6768E  02 

5.002  7E-01 

6.08756-01 

1.-40 

4 . 7439E-0 l 

2.6398E  02 

4 . 8460E-01 

5.6906E-01 

1.50 

4.7434E-C1 

2.6241E  C2 

4.8708E-01 

5. 5656E-01 

1.60 

4.5536E-01 

2.5948E  02 

4.9687E-01 

5.4500E-01 

1.70 

4.3167E-01 

2. 57766  02 

5.1249E-01 

5.4406E-01 

l  .«0 

4.Q524E-01 

2.5254E  02 

5.1803C-01 

4. 7531E-01 

1.90 

3.9663E-01 

2.5162E  02 

5.2362E-01 

4 . 74C6E-01 

2.  0 

2. 1C 

** 

2.20 

** 

2.30 

** 

2.40 

** 

2.50 

** 

2.60 

** 

2.70 

** 

2.80 

** 

2.90 

4.02 1 6E-0 1 

2.6078E  02 

5.305UE-01 

5.4406E-01 

NC  S'.AP-THROUGH 


45 


AFFDL-TR-68-136 


TABLE  10. 


BUCKLING  LOADS  FOR  IMPERFECT  GLASS-LPOXY 
COMPOSITE  CYLINDRICAL  SHELL  FOR 

0.30 


CASE 

<r 

max 

ar 

max 

max 

I.  0 

** 

1.10 

*♦ 

1.20 

** 

1.30 

** 

1.40 

** 

1.50 

** 

1.60 

** 

1.70 

♦* 

1.80 

3. 3696E-01 

2.0999E 

02 

4.3075E-01 

4.6313E-01 

1.90 

3.3176E-C1 

2.1047E 

02 

4.3798E-01 

4.7531E-01 

2 .  0  ♦♦ 

2.10  ** 

2.20  ** 

2.30  ** 

2.4 0  ♦* 

2.50  #* 

2.60  ** 

2.70  ♦* 

2.80  ♦* 

2.90  *♦ 

**  NC  SNAP- ThRCUGH 


46 


AFFDL-TR-68-136 


TABLE  11. 

BUCKLING  LOADS  FOR  IMPERFECT  BORON-EPOXY 
COMPOSITE  CYLINDRICAL  SHELL  FOR 

O.Ol 


CASE 

amax 

nr 

max 

P* 

£mox 

1.  0 

4.65426-01 

7.8144E 

02 

9.4230E-01 

4. 76256-01 

1.10 

5. 22876-01 

8.5677E 

02 

9.2996E— 01 

5. 3187E-01 

1.20 

6.6635E-01 

1.0142E 

03 

9 . 1092E-01 

6.73176-01 

1.30 

8.47896-01 

1.1660E 

03 

8.99846-01 

8.5292E-01 

1.40 

9. 2045E-01 

1.1960E 

03 

9.0524E— 01 

9.2442E— 01 

1.50 

8.3492E-01 

L.1097E 

03 

9. 1233E-01 

8.39926-01 

1.60 

6.9680E— 01  * 

9.97716 

02 

9.1676E-01 

4.4200E-03 

1.70 

5.6065E-01 

8.9306E 

02 

9. 28636-01 

5.7531E-01 

*1.80 

4. 3423E-01 

7.98816 

02 

9.4672E— 01 

4.6000E-01 

'  1.90 

T.63576-01 

7.43C9E 

02 

9.5992E-01 

4.0438E-01 

2.  0 

3. 5321E-01 

7.4827E 

02 

9.6661E-01 

4.4156E-01 

2.10 

3.8825E-01 

7.9572E 

02 

9.4304E-01 

4. 1469E-01 

?.20 

4.7535E-01 

8.8240E 

02 

9. 1754E-01 

4.91876-01 

2.30 

6.2007E-01 

9.8811E 

02 

9.0794E-01 

6.37816-01 

2.40 

8.11356-01 

1.11576 

03 

9.17206-01 

8. 35946-01 

2.50 

9.31246-01 

l .21216 

03 

9*17396-01 

9.5125E-01 

2.60 

8.65546-01 

1.1746E 

03 

9.0646E-01 

8.77666-01 

2, 70 

6.93386-01 _ 

1.01 74 E 

03 

9, 1383E-01 

7.01346-01 

2. 80 

5.3844E-01 

8.58166 

02 

9.31466-01 

5.48126-01 

2.90 

4.6542E-01 

7.81446 

02 

9. 42306-01 

4.76256-01 

AFFDL-TR-68-136 


TABLE  12. 

BUCKLING  LOADS  FOR  IMPERFECT  BORON-EPOXY 
COMPOSITE  CYLINDRICAL  SHELL  FOR 

¥*■=  0.02 


CASE 

°ma« 

a  max 

/>* 

1.  0 

4.4868E-01 

7.5333E 

02 

9.0841E-01 

4. 67506-01 

1.10 

5.C030E-01 

8.1979E 

02 

8 . 8982E-0 1 

5.1656E-01 

1.20 

6. 3G49E-01 

9.5964E 

02 

8.6190E— 01 

6.42506-01 

1.30 

7.9694E-01 

1.0959E 

03 

8.4577E-01 

6.4200E-03 

1.40 

8.6638E-01 

1.1257E 

03 

8.5207E-01 

7.3670E-02 

1.50 

7 . 8865E-0 1 

1.0482E 

03 

8.6177E-01 

7.9717E-01 

1.60 

6 . 59256-01 

9.4394E 

02 

8.6735E-01 

6.69846-01 

1.70 

5.35756-01 

8.5340E 

02 

8.8739E-01 

5. 61566-01 

1.80 

4.2029E-01 

7.73L7E 

02 

9 . 1632E-0 1 

4.68756-01 

1.90 

3. 55056-01 

7.2568E 

02 

9.3743E-01 

4.35316-01 

2.  0 

3.4655E-01 

7.3417E 

02 

9.4839E-0i 

4.89C6E-01 

2.10 

3.7532E-01 

7.6922E 

02 

9 . 1 163E-01 

4.28446-01 

2.20 

4.5249E-01 

8.3996E 

02 

8. 73416-01 

4.85006-01 

2.30 

5 .87906-0 1 

9.3685E 

02 

8 .60846-01 

6.24C6E-01 

2.40 

7.7401E-01 

1.06436 

03 

8.7499E-01 

8.25636-01 

2.50 

8.8743E-01 

1.15506 

03 

8. 74236-01 

9. 27 19E-01 

2.60 

8. 1911E-01 

1.1116E 

03 

8 • 5783E-01 

8.39376-01 

2.70 

6.5743E-01 

9.6462E 

02 

8 . 6645E-01 

6. 72196-01 

2.80 

5. 1568E-01 

8.2189E 

02 

8.9209E-01 

5.32816-01 

2.90 

4.4868E-01 

7.5333E 

02 

9.0841E-01 

4.67506-01 
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TABLE  13  e 

BUCKLING  LOAOS  FOR  IMPERFECT  BORON-EPOXY 
COMPOSITE  CYLINDRICAL  SHELL  FOR 

VTS=  0.04 


1 


CASE 

&  max 

pm 

B  max 

1.  0 

4.2315E-01 

7.1047E 

02 

8 . 5672E-01 

4. 5719E-01 

1.10 

4.66736-01 

7.6478E 

02 

8.301 IE-01 

4. 97196-01 

1.20 

5.7881E-01 

8 . 8058E 

02 

7.9125E-01 

6.02 19E-01 

1.30 

7.2419E-01 

9.95906 

02 

7.6856E-01 

7.40766-01 

1.40 

7 .67306-01 

1 .02306 

03 

7.7429F-01 

0.0092E-01 

1.50 

7.1953E-01 

9.56316 

02 

7.86246-01 

7.3542E-01 

1 .60 

6.0485E-01 

8.6605E 

02 

Y.9578E-01 

6.2844E-01 

1.70 

4.9885E— 01 

7.94626 

02 

8.2627E-01 

5.48756-01 

1.80 

4.0034E-01 

7.3647E 

02 

8. 7283E-01 

5. 1531E-01 

1.90 

3.42276-01 

6.9956E 

02 

9.0368E-01 

4.9531E-01 

2.  0 

3. 3537E-01 

7.1048E 

02 

9. 1779E-01 

5.428IE-01 

2.10 

3.5769E-01 

7.33C9E 

02 

8.6881E-01 

5.0562E-01 

2.20 

4.2141E-01 

7.8227E 

02 

8. 1342E-01 

5. 00626-01 

2.30 

5.4551E-01 

8.6930E 

02 

7.9877E-01 

6.38446-01 

2.40 

** 

2.50 

8.2836E-01 

1 .07816 

03 

8. 1604E-01 

9.24066-01 

2.60 

7 . 5430E-0 1 

1.0236E 

03 

7 . 8996E-01 

7.9687E-01 

2.70 

6.0571E-01 

8.8874E 

02 

7.9829E-01 

6*  3469E-01 

2.80 

4^81786-01 

7 .67866 

02 

8.3344E-01 

5. 1406E-01 

2.90 

4.2315E-01 

7.10476 

02 

8.56726-01 

4.5  7 19E-01 

**  NO  SNAP-THROUGH 
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TABLE  14. 


8UCKLING  LOADS  FOR  IMPERFECT  BGKCN-SPOXV 
COMPOSITE  CYLINDRICAL  SHELL  FOR 

w“  =  0.06 

CASE 

°mar 

P* 

£max 

1.  0 

4.0299E-01 

6.7662E  02 

8. 1590E-01 

4.52816-01 

i.io 

4.4093E-01 

7.2251E  02 

7.8422E-01 

4.8656E -01 

1.20 

5.4018E-01 

8 .22186  02 

7.38456-01 

5.76256-01 

1.30 

6.7089E-01 

9.22606  02 

7. 1199E-0I 

6.99O6E-01 

1.40 

7.1765E-01 

9.3249E  02 

7.0579E-0 1 

7.5098E-01 

1.50 

6.6619E-01 

8.85416  02 

7.2796E-01 

6.9C31E-01 

1.60 

5.6222E-01 

8. 050 IE  02 

7 • 3970E-01 

5.9687E-01 

1.70 

4,71136-01 

7 .50466  02 

7.8035E-01 

5.5187E-01 

1.80 

#*- 

1.00 

3.3184E-0L 

6.7824E  02 

8.7615E-01 

5.4281E-01 

2.  0 

3.2564E-01 

6.09876  02 

8.9116E-01 

5.8760E-01 

2.10 

** 

2.20 

2.30 

** 

2.40 

>»« 

2.50 

** 

2.60 

7.0797E-01 

9.6076E  02 

T.4144E-01 

7. 78  75E-01 

2.70 

5.6720E-01 

8.32236  02 

7.4754E-01 

6.12196-01 

2.00 

4.5570E-01 

7,26296  02 

7.8833E-01 

5.C437F-01 

2.90 

4.0299E-01 

6.76626  02 

8.15906-01 

4.5281E-01 

N C  SNAP- 

THROUGH 
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TABLE  15. 

BUCKLING  LOADS  FOR  IMPERFECT  BORON-EPOXY 
COMPOSITE  CYLINDRICAL  SHELL  FOR 

0.08 


CASE 

^mos 

/>* 

1,  0 

3. 86026-01 

6.48136  02 

7. 81546-01 

4^52816-01 

1.10 

4. 1975E-01 

6.S780E  02 

? . 4655E-01 

4.8281E-01 

1.20 

5.09S6E-0I 

7.7482E  02 

6 . 9590E-01 

5.5937E-01 

1.30 

6. 2804E-0 1 

8.6367E  02 

6.6652E-01 

6.6719E-0I 

1.40 

6.78820-01 

8.Q204E  02 

6.67606-01 

7.0594E-01 

1.50 

6. 22526-01 

8.2737 E  02 

6.8024E-01 

6. 53126-01 

X.6C 

5.27211-01 

7.5488E  02 

6.9363E-0I 

5.7406E-01 

lo?0 

4.4974E— 0i~ 

7.16396  02 

7.4492E-01 

5.8844E-01 

1.80 

$4 

1.90 

3.2270E-0I 

6o5956E  02 

8.5201E— 01 

5.8644E-01 

2.  0 

3. 1696E-01 

6.7148E  02 

0.6741E-O1 

6.35006-01 

2.10 

2.20 

2.30 

2.40 

2.50 

2.60 

6.72T0E-01 

9.1289E  0;j 

7.0450E-01 

7.85946-01 

©• 

» 

rv 

5 » 3&36E-01 

7.8698E  02 

7*06896-01 

6.0000E-01 

2.80 

4.34276-01 

6. 92 146  02 

7.5125E-01 

5.0156E-01 

2.90 

3.8602E-01 

6.4813E  02 

7.81546-01 

4.5281E-01 

**  NO  SNAP-THROUGH 
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TABLE  16. 


BUCKLING  LOADS  FOR  IMPERFECT  BORON-EPOXY 
COMPOSITE  CYLINDRICAL  SHELL  FOR 

is=  0.10 


CASE 

^mca 

&  mem 

P 

^mox 

l.  0 

3. T129E-01 

6.2339E 

02 

7.5172E-01 

4.5625E-01 

1.10 

4. 01795-01 

6.5837E 

02 

7 . 14615-01 

4.8469E— 01 

1.20 

4.83I0E-01 

7.3531E 

02 

6.6041E-01 

5.4969E— 01 

1.30 

5.9248E-01 

8.1477E 

02 

6.2878E-01 

6.4406E-G1 

1.40 

6.3786E-01 

8.2S82E 

02 

6.2732E-01 

6.7281E-01 

1.50 

5.3541E-01 

7.7805E 

02 

6.3969E— 01 

6.2437E-01 

1 . 60 

4.9766E-01 

7.1257E 

02 

6.5476E— 01 

5.5844E-01 

1.70 

** 

I. 30 

** 

1*90 

3.1454E-01 

6.4288E 

02 

8.3047E-01 

6.3875E-01 

2.  0 

3.0918E-01 

6.5500E 

02 

8.4612E-01 

6.9250E-01 

2.10 

2.20 

o 

0 

(V 

** 

2.40 

2.50 

2.60 

** 

2.70 

5. 1 985E-01 

7.49555 

02 

6.7327E-0I 

5.9719E-01 

2.30 

4.1614E-01 

5.6324E 

02 

7 . 1989E-01 

5.05316-01 

2.90 

3.7129fe-01 

6.2339E 

02 

7.5172E-01 

4. 56256-01 

**  NO  SNAP-THROUGH 
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TABLE  17. 

BUCKLING  LOADS  FOR  IMPERFECT  BORON-EPOXY 
COMPOSITE  CYLINDRICAL  ShELL  FOR 

0.20 


CASE 

^ma* 

tfmaii 

£max 

1.  0 

3. 1791E-C1 

5.3377E  02 

6.4365E-01 

5.2500E-01 

l  .10 

1.20 

** 

1.30 

** 

1.40 

5.02 12E-01 

6.5244E  02 

4.9382L-01 

6.0ie7E-01 

1.50 

4.6213E— 01 

6.1420E  02 

5.0498E-01 

5.8187F-01 

1.60 

**> 

1.70 

** 

1.80 

** 

P-* 

0 

4D 

O 

** 

2.  0 

2.10 

<■* 

2,20 

2.30 

♦  d 

2.40 

** 

2.50 

2.60 

** 

2 . 70 

2.80 

. 

2.90 

3.1791E-01 

5.3377E  02 

6.4365E-01 

5.2500E-01 

**  NO  SNAP-THROUGH 
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Nl3  ABSTRACT 
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The  effect  of  initial  geometric  imperfections  on  the  buckling  and  poetbuckling  behavior 
of  composite  cylindrical  shells  subjected  to  uniform  axial  compression  is  studied  in  this 
report.  The  solution  is  obtained  by  employing  von  Karm&n-Donnell  nonlinear  strain- 
displacement  relations  and  the  principle  of  stationary  potential  energy.  Numerical  results 
are  given  for  various  fiber  orientations  in  the  three-layer  shell  consisting  of  either 
glass-epoxy  or  boron-epoxy  composites,  with  different  initial  imperfections.  Results 
indicate  that  the  boron-epoxy  composite  shells  are  less  imperfection  sensitive  than  the 
glass -epoxy  composite  shells.  Isotropic  shells  are  found  to  be  more  imperfection  sensitive 
than  composite  shells.  It  is  noticed  that  the  increase  or  decrease  in  the  classical  buckling 
load  with  change  in  fiber  orientation  is  generally  accompanied  by  a  decrease  or  increase  in 
imperfection  sensitivity  of  the  shell. 
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